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The recently obtained complete solution of the simultaneous diagonalization of matrices H A and H in the 
hydrodynamic diffusion equation has basically changed the diagonal values v v of the symmetric matrix H of 
hydrodynamic interaction between all the heads of the elastic random coil model of the isolated macromolecule in 
solution. Since these values enter explicitly the expressions for the intrinsic stress and refractive index tensor in an 
alternating flow field if based on the concept of internal viscosity of the model one had to recalculate all values 
obtained formerly by using the then generally accepted erroneous set of v v data. The new v p equal unity 
independent ofp while the old values were larger than 1 for small/) and smaller for large p. Hence their too large 
contributions in the former range are partially compensated by their too small contributions in the latter region. As 
a consequence in the whole range investigated, between 3 and 300 chain links, the differences in theological and 
rheooptical effects are relatively small, up to a factor of 2, although at higher link number the differences tend to 
grow with the logarithm of this number. 
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1. Introduction 

The correct simultaneous diagonalization [1J 1 of H A and 
H matrices in Zimm's hydrodynamic equation [2] for the 
ideally flexible necklace model of randomly coiled isolated 
linear macromolecule in laminar flow makes possible a more 
adequate calculation of intrinsic stress tensor in all those 
cases where the coil is not yet noticeably deformed by the 
flow. Such a zero gradient case includes the frequency de- 
pendence of viscosity [2], shear modulus, shear birefrin- 
gence, normal stress difference, and acoustic birefringence 
but not the gradient dependence of these effects. 

The main change introduced by the new solution as com- 
pared witb the older incomplete solutions [3-8] is not in the 
eigenvalues k p of H A 



Q 1 HAQ 



(1) 



which were already calculated correctly in recent papers 
[6-8] and even tabulated for Z between 1 and 15 and h* = 
(3/tt) 1/2 a h /b = .01, .1, .2, .3, [6] and for Z = 250, h* = 
.3 andZ = 300, h* = 0.4 [1]. Here ah is the hydrodynamic 
radius of the bead, Z + 1 is the number of beads, Z is the 
number and bo the root mean square length of the links. The 



difference is in the diagonal elements v p of the tensor 

Q 1 H Q 17 = N = 1 

which turns out to be the unity tensor with all v p — 1 . 
makes the diagonal elements of 

M = Q 7 AQ = A/N = A 



(2) 
Phis 

(3) 



agree with those of A, i.e. \x p = \ p . Here Q is the transfor- 
mation matrix of the original 3(Z +1) dimensional vector r 
of bead coordinates to that of dimensionless normal coordi- 
nates u 



r = b Q u 



(4) 



and Q 7 its transpose. 

Equation 3 completely differs from the original estimate [3] 
that in first approximation the diagonal elements of the matrix 
M equal the eigenvalues \ p r of the Rouse model with vanish- 
ing hydrodynamic interaction, i.e. 






4 sin 2 [pir/2(Z + 1 )] 



and hence 



1 Figures in brackets indicate the literature references at the end of this paper. 



VpZ — ^pZ/^pR- 



(5) 



(6) 
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Figure 1. Diagonal elements v p (new theory) and Pp } = 
theory, broken lines) for Z = 100 and h* = 0, .1, .2, 

For h* = (Rouse) the values of both theories coincide. 



100 



^•pzMpr (old 
.3, .4. 



The subscripts Z and R refer to Zimm and Rouse [9] model, 
respectively. Such an estimate was based [3] on the supposi- 
tion that the transformation matrix Q changes so little by the 
introduction of hydrodynamic interaction that M remains 
practically the same as in the free draining case, i.e. M = 

A* 

The diagonal elements v p io\Z = 100 and h* = (Rouse), 
.1, .2, .3, and .4 (Zimm) are plotted in figure 1 and the 
eigenvalues k p collected in table 1. The old values of v p l) — 
kpz/^pR are connected with a broken line. They are partially 
situated below and partially above v p = 1. This reduces the 
differences between the old and new values in quantities 
dependent on v p . They show up in the excess stress tensor as 
soon as internal viscosity, defined by the frictional parameter 
<£, is explicitly considered. The very large old values v p for 
small p do not matter very much because in all formulae they 
are multiplied by very small values (f p = p(p/Z . Since the 
differences between old and new v p increase with /i*, in that 
which follows, the comparison of calculated effects will be 
mainly done for/i* = .4, i.e., for a very large hydrodynamic 
interaction. The differences are smaller for smaller h* and of 
course disappear for the free draining coil with h* = where 
kpz — kpR- 

The new values /jl p and v p do not enter the conventional 
expressions for the intrinsic stress or birefringence tensor of 
perfectly flexible necklace model so that no changes occur in 
these quantities. The situation is completely different if one 
considers the effects of internal viscosity which depend on v p 
[10-26]. They will be most conspicuous in the values of 
viscosity corresponding to high frequency flow field and in 
the phase angle between the stress and strain rate or between 
the birefringence and strain rate. 

The interaction tensor H depends on the inverse intrabead 
distances l/r^ which makes the hydrodynamic diffusion 
equation intrinsically non-linear. By replacing l/r jk with its 
average value (l/r^) the tensor H becomes a constant which 
makes the diffusion equation linear and hence allows the 
introduction of normal coordinates according to eqs (1) to (4). 



Such a procedure eliminates the possibility of any realistic 
consideration of the gradient dependence of any rheological 
or rheooptical effect because it does not take into account the 
change of shape of the random coil in flow which expands the 
molecule and hence increases the interbead distances r jk . 
Note also that by preaveraging over all angles between the 
velocity and the interbead vector the formulation of H as 
function of l/r^ completely evades the consideration of ani- 
sotropy of hydrodynamic interaction which by itself yields a 
gradient dependence of intrinsic viscosity [27] large enough 
for explaining experimental data. 

The general toleration of such a profound modification of 
hydrodynamic interaction by the replacement of I /r jk with its 
average (l/r jk ) makes hard to understand the almost general 
objections to the introduction of internal viscosity as a resist- 
ance of the necklace model against the deformational compo- 
nent of the normal eigenmodes [28-30]. If one accepts the 
rather questionable linearization of the hydrodynamic diffu- 
sion equation one has to accept also the next step, i.e., the 
concept of internal viscosity based on this linearity and its 
introduction in such a manner that the mathematical treat- 
ment remains as simple as possible. 

In that which follows the results of the new theory will be 
compared with those of the old one for h* = .4 and Z — 100 
in the whole frequency range and the dependence of the 
limiting values for co = cc on Z in the range between Z — 3 
and 300 and on /i* in the range between .1 and .4. In all 
cases the ratio between the internal viscosity coefficient <p 
and the frictional coefficient of the bead/ = 67Ta n r) s will be 
assumed constant, <p/f = 2. Here r) s is viscosity of the 
solvent. The subscript s applies to the properties of the 
solvent. The corresponding non-subscripted quantities relate 
to solution. 

2. Internal Viscosity 

The concept of internal viscosity was introduced in order to 
express the inability of the randomly coiled polymer molecule 
to change rapidly its shape [31, 32, 10, 11]. Such changes 
occur during the rotation of the macromolecule in a flow with 
a rotational component, e.g., the laminar flow with transverse 
gradient, when the individual segments are alternatively 
passing from the direction of compression to the direction of 
extension and vice versa. The directions of maximum 
compression and extension of the volume element are in the 
flow plane perpendicular to each other. The rapidity of 
change is given by the transverse gradient which equals twice 
the angular velocity of the ideally flexible coil which rotates 
with the volume element. The other case is the oscillating 
flow field where the oscillation frequency determines the 
rapidity of change from compression to extension and vice 
versa. 

In the limiting cases of zero gradient and zero frequency 
the deformation inability of the macromolecule does not play 
any role. All the changes occur so slowly that the effects are 
identical for completely rigid and ideally flexible coils if only 
their conformational distributions agree with each other. 
With increasing gradient and/or frequency, however, the 
time effects are playing a gradually increasing role. They are 
maximum in the second Newtonian regime corresponding to J 
—> °° or CO — > oo. 

The rigidity of the macromolecule can be assigned to 
different properties of the chain. It can be caused by the 
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Table 1. Eigenvalues A. 


p /orZ = 100 and h* = 


(Rouse), 


.7, .2, .3, A (Zimm) 






p 


//* = 


.1 


.2 


.3 


.4 


P 


h* = 


.1 


.2 


.3 


.4 


1 


0.000967 


0.002170 


0.003345 


0.004506 


0.005643 


51 


2.03111 


1.7822 


1.5334 


1.2846 


1.0359 


2 


.003869 


.007318 


.010718 


.014084 


.017417 


52 


2.09329 


1.8320 


1.5707 


1.3094 


1.0483 


3 


.008701 


.014650 


.020541 


.026383 


.032173 


53 


2.15538 


1.8815 


1.6077 


1.3339 


L.0602 


4 


.015460 


.023986 


.032447 


.040848 


.049177 


54 


2.21731 


1.9308 


1.6443 


1.3579 


1.07k) 


5 


.024139 


.035175 


.046144 


.057043 


.067854 


55 


2.27904 


1.9799 


1.6807 


1.3816 


1.0826 


6 


.034730 


.048140 


.061482 


.074746 


.087907 


56 


2.34050 


2.0286 


1.7167 


1.4049 


1.0932 


7 


.047220 


.062802 


.078315 


.093745 


.10906 


57 


2.40153 


2.0770 


1.7523 


1.4277 


1.1033 


8 


.061602 


.079111 


.096552 


.11391 


.131113 


58 


2.46236 


2.1249 


1.7875 


1.4502 


1.1130 


9 


.077859 


.097010 


.11610 


. 13509 


. 15395 


59 


2.52265 


2.1725 


1.8223 


1.4722 


1.1223 


10 


.095974 


.11646 


.13688 


M5721 


.17740 


60 


2.58244 


2.2195 


1.8566 


1.4938 


1.1311 


11 


.115932 


.13742 


.15884 


.18017 


.20135 


61 


2.64166 


2.2661 


1.8905 


1.5150 


1.1396 


12 


.137713 


.15984 


.18192 


.20390 


.22571 


62 


2.70026 


2.3120 


1.9238 


1.5357 


1.1477 


13 


.161295 


.18370 


.20605 


.22831 


.25040 


63 


2.75819 


2.3574 


1.9566 


1.5559 


1.1554 


14 


. 186656 


.20895 


.23120 


.25335 


.27533 


64 


2.81538 


2.4021 


1.9889 


1.5757 


1.1627 


15 


.213771 


.23556 


.25730 


.27896 


.30044 


65 


2.87178 


2.4462 


2.0206 


1.5950 


1.1697 


16 


.242614 


.26350 


.28434 


.30509 


.32566 


66 


2.92734 


2.4895 


2.0517 


1.6139 


1 . 1 763 


17 


.273158 


.29272 


.31224 


.33169 


.35095 


67 


2.98200 


2.5321 


2.0821 


1.6323 


1.1826 


18 


.305372 


.32320 


.31099 


.35870 


.37()23 


68 


3.03571 


2.5738 


2.1120 


1.6502 


1.1885 


19 


.339225 


.35490 


.37053 


.38010 


. 10148 


69 


3.08842 


2.0148 


2.1411 


1.6676 


1.1941 


20 


.374681 


.38778 


.40084 


.41383 


.42664 


70 


3.14007 


2.6548 


2. 100(> 


1.6845 


1.1995 


21 


.411719 


.42180 


.43186 


.44186 


.45168 


71 


3.19063 


2.6940 


2.1974 


1.7009 


1.2045 


22 


.450288 


.45694 


.46358 


.47015 


.47655 


72 


3.24003 


2.7323 


2.2245 


1.7168 


1.2092 


23 


.490356 


.49316 


.49594 


.49867 


.50123 


73 


3.28823 


2.7695 


2.2508 


1.7322 


1.2137 


24 


.531886 


.53041 


.52892 


.52738 


.52568 


74 


3.33518 


2.8058 


2.2764 


1.7471 


L.2179 


25 


.574835 


.56867 


.56248 


.55626 


.54987 


75 


3.38085 


2.8410 


2.3012 


1.7614 


1.2218 


26 


.619163 


.60788 


.59659 


.58527 


.57378 


76 


3.42518 


2.8752 


2.3252 


1.7753 


1.2255 


27 


.664827 


.64802 


.63122 


.61438 


.59738 


77 


3.46812 


2.9083 


2.3484 


1.7886 


1.2290 


28 


.711782 


.68905 


.66632 


.64356 


.62065 


78 


3.50965 


2.9402 


2.3708 


1.801 1 


1.2322 


29 


.759984 


.73093 


.70187 


.0727') 


.(.1357 


79 


3.54972 


2.9710 


2.3924 


1.8137 


1.2352 


30 


.809386 


.77361 


.73784 


.70205 


.66613 


80 


3.58829 


3.0007 


2.4131 


1.8255 


1.2380 


31 


.859939 


.81706 


.77419 


.73130 


.68829 


81 


3.62532 


3.0291 


2. 1329 


1.8367 


1.2406 


32 


.911596 


.86124 


.81089 


.76053 


.71005 


82 


3.66078 


3.0563 


2.4518 


1.8474 


1.2430 


33 


.964305 


.90610 


.84790 


.78970 


.73140 


83 


3.69464 


3.0823 


2.4699 


1.8575 


1.2453 


34 


1.01802 


.95160 


.88521 


.81881 


.75231 


84 


3.72685 


3.1069 


2.4870 


1.8672 


1.2474 


35 


1.07268 


.99771 


.92277 


.84782 


.77279 


85 


3.75737 


3.1303 


2.5033 


1.8762 


1.2493 


36 


1.12824 


1.0444 


.96055 


.87672 


.79281 


86 


3.78624 


3.1524 


2.5186 


1.8848 


1.2510 


37 


1 . 18464 


1.0916 


.99852 


.90548 


.81238 


87 


3.81335 


3.1731 


2.532<) 


1.8928 


1.2526 


38 


1.24183 


1.1392 


1.0367 


.93410 


.83148 


88 


3.83871 


3.1925 


2.5464 


1.9002 


1.2541 


39 


1.29975 


1.1873 


1.0749 


.96254 


.85011 


89 


3.86229 


3.2105 


2.5588 


1.9071 


1.2554 


40 


1.35835 


1.2358 


1.1133 


.99079 


.86826 


90 


3.88407 


3.2272 


2.5703 


1.9134 


1.2566 


41 


1.41758 


1.2846 


1.1517 


1.0188 


.88592 


91 


3.90403 


3.2424 


2.5808 


1.9192 


1.2577 


42 


1.47736 


1.3338 


1.1902 


1.0467 


.90311 


92 


3.92215 


3.2563 


2.5904 


1.9245 


1.2586 


43 


1.53758 


1.3832 


1.2287 


1.0742 


.91980 


93 


3.03840 


3.2687 


2.5989 


1.9292 


1.2595 


44 


1.59839 


1.4328 


1.2671 


1.1016 


.93601 


94 


3.95278 


3.2796 


2.6065 


1.9333 


1.2602 


45 


1.65952 


1.4825 


1.3056 


1.1286 


.95173 


95 


3.96527 


3.2892 


2.6130 


1.9369 


1.2608 


46 


1.72098 


1.5324 


1.3439 


1.1554 


.96696 


96 


^.97586 


3.2972 


2.6186 


1.9400 


1.2614 


47 


1.78270 


1.5824 


1.3821 


1.1819 


.98171 


97 


3.98454 


3.3038 


2.6232 


1.9425 


1.2618 


48 


1.84464 


1.6324 


1.4202 


1.2021 


.99597 


98 


3.99130 


3.3090 


2.6267 


1.9444 


1.2621 


49 


1.90673 


1.6824 


1.4581 


1.2339 


1.0098 


99 


3.99613 


3.3127 


2.6292 


1.9458 


1.2624 


50 


1.96890 


1.7324 


1.4959 


1.2594 


1.0231 


100 


3.99903 


3.3149 


2.6308 


1.9466 


1.2625 



energy barrier separating the gauche and trans conformations 
which makes any conformational change more time consum- 
ing than in the case of perfectly soft model [10]. Since the 
height of the energy barrier is independent of the viscosity of 
the solvent its relative effect as measured by the ratio of 
internal frictional resistance (p caused by the barrier and the 
frictional resistance / of the segment decreases with the 
viscosity of the solvent. The macromolecule acts as very rigid 
in a low viscosity liquid, e.g., acetone with 7) s — .322cPoise 
= 32.2 mNs/m 3 , and as very flexible in a high viscosity 
solvent, e.g., Aroclor with 7] s up to 100 Poise (= 1 kNs/m 3 ) 
and higher. 



Another cause of slow molecular response to the rapidly 
changing flow field resides with the conformational restraints 
of the chain which permit only an interchange of gauche and 
trans conformations [20, 23, 24]. With almost rigid length of 
valency bonds this means that most changes of length and 
position of any chain segment require a much larger segment 
displacement than formulated in the ideally flexible necklace 
model which does not consider any inherent limitation of 
bead motion. Generally an axial displacement of the segment 
requires also some lateral displacement and vice versa. As a 
consequence the resistance of beads to position change is 
larger than assumed on the basis of hydrodynamic radius a^ 
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The ratio of the so obtained coefficient <p to/ is independent 
of solvent viscosity because (p and/ are both proportional to 
r) s . Their ratio just measures the ratio of true displacement to 
the minimum displacement explicitely considered in the 
diffusion equation. It seems to be close to 2 for vinyl poly- 
mers 

The effect of internal viscosity is formulated in the system 
of normal modes [10]. For the pth normal mode of deforma- 
tion one has a resistance coefficient (pp/Z. Such a choice is 
reasonable for both origins of internal viscosity as just dis- 
cussed. In the first case one can argue that the changes to 
comply with any deformational mode are linearly increasing 
with the number of chain atoms between subsequent nodes, 
i.e., with Z/p, because a conformational change takes place 
with equal probability at any of these chain atoms. This 
makes the resistance increase with p/Z. In the second case 
the displacement at lower modes can be achieved in many 
ways so that the actual lengthening of displacement path is 
much less noticeable than at higher modes where the confor- 
mational restrictions are soon becoming of overwhelming 
importance. 

One may argue that the whole concept of internal viscosity 
can be discarded because it is not based on some strictly 
fundamental analysis of chain dynamics. It was indeed intro- 
duced in a rather pragmatic manner which also permitted an 
easy mathematical treatment [10]. But it turns out that all 
more detailed treatments of Brownian motion of beads or of 
correlation between the motion of two or more beads [33-35] 
lead to some, often hidden, statement of molecular rigidity 
which is needed for the results of such a study to reproduce 
the characteristic rheological features of polymer systems, 
e.g., the non-vanishing limiting intrinsic viscosity at very 
high frequency [36-51]. Such a state of affair seems more to 
support than to refute the concept of internal viscosity in spite 
of its more pragmatic than fundamental way of introduction. 

3. The Distribution Function of the Beads 

The continuity equation of the ideally flexible necklace in 
laminar flow which determines the distribution function i//(r , 
r i * * * r z) reads 

V/ [vijj - (3D /b 2 ) H A r i// - D H V r ifl 

= -difj/dt. (7) 

Here D = kT/f is the translational diffusion coefficient of 
the bead. By introduction of normal coordinates, eq (4), one 
transforms eq (7) into a system of Z + 1 partial differential 
equations 

V P [v p 4/ p ~ (SD /b 2 )k p u p - (D /b 2 )v p V p \jj p ] 

= -dtyjdt (8) 

each depending only on the coordinates of the/?th eigenmode. 
Note that Vp in normal coordinates has the dimension s -1 . 
The distribution function of the coil is the product of all ifj p 



i//(u, t) = tWuj, t) • • ■ i// z (u z , t). 



(9) 



The functions \fj p depend on the kind of flow field v. The 0th 
mode does not show up in i// because it represents a uniform 
translation of the whole necklace which does not affect i//. 



The introduction of internal viscosity adds a viscous type 
resistance coefficient <p p = p<p/Z opposing thepth eigenmode 
of the true deformation rate of the coil. This rate is obtained 
by subraction of pure rotational velocity Cl X u p from the 
total deformation rate du p /dt. This yields an internal viscos- 
ity force [10, 11] 

Fip = -(p<p/Z)(du p /dt -ax u p )b . (10) 

If one introduces this force in the pth diffusion equation (eq 
(8)) one obtains after some rearrangements 

(i + vtfp/fr 1 v P {[v p - (3zv&o 2 ) v* + (v P <pp/j)n 

x Up - (zVVVpVpMp) = -dxjj/dt. (li) 

The distribution function i// p depends on the kind of flow and 
on the angular velocity vertor O. 

In a jet or plane flow with longitudinal gradient without a 
rotational flow component one has Q = 0. In a flow with 
transverse gradient, v = y(y,o,o), the angular velocity 
equals —yl'2 



ft = (0, 0, - y/2) 



(12) 



for relatively soft molecules which rotate in phase with the 
volume element. This is the case with practically all conven- 
tional macromolecules if the degree of polymerization is so 
high that a truly random coil is formed. Very short chains, 
chains with a great many double bonds, ladder type and 
multiple strand molecules, however, are more rigid and tend 
to rotate with a non-uniform angular velocity which depends 
on orientation of the molecule. It is different from that of he 
volume element. Assymptotically, at very high rigidity and 
fully extended shape of the macromolecule, it approaches 
that of rigid bodies, e.g., rods or ellipsoides. In that which 
follows only the case of practically undeformed relatively soft 
coils with fl = — yl'2 will be considered. 

In the steady state flow with transverse gradient the p\h 
eigenmode distribution function of the soft necklace reads 



*M£p, Vp, €p) 



exp 



IT 



i + & 2 -/V y /2 



i+/$ 



^{(1 + iS; 2 - ft, 2 )& 2 - 2j3 p &t ?p 



+ [i + p'p (&> + ft,)}*?, 2 ) +■ (i + P' p 2 Kp 2 } 

Pp = yb 2 /6D k p 
P' p = j8 p (l + v p <p p /f). (13) 



The index p runs from 1 to Z . The value is excluded. 
Without internal viscosity, <p p — 0, one has /3' p = /3 P and one 
obtains the conventional distribution function. Note that eq 
(13) and the distribution function eq (6.1) in Ref. [1] refer to 
different flow fields: 7(7,0,0) in the former and 7(2,0,0) in the 
latter case. 
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Equation (13) explicitely contains the diagonal terms \ p , 
}JL P and v p . In the averages (dyadic formulation) 

z 
<r T Ar) = V fc T Mu ) = 6 2 2 iM"^* > (14) 



which appear in the intrinsic stress and optical tensor the 
coefficients /x p drop completely because the averages (upUp) 
are proportional to l//x p . The coefficients v p , however, re- 
main in /3' p as long as <p =£ 0. As already mentioned they are 
equal to unity (new theory) while formerly they were approxi- 
mated by (AzAtf)p = k P (h*)/k p (h* = 0) ± 1 (old theory). 

In the oscillating flow field the gradient is a function of 
time 



7o^ 



(15) 



The amplitude y is so small that the molecule remains 
practically undeformed so that the zero gradient eigenvalues 
\ p and diagonal elements v p , calculated for the coil at rest, 
are applicable. 

4. The Intrinsic Stress Tensor 

In dyadic formulation the intrinsic stress tensor reads 



Here N is Avogadro number, M is molecular weight, and F is 
the vector of forces exerted by the beads of the necklace 
model on the flowing liquid. One has in the space of normal 
coordinates 

F = (1 + N^/Z^Cr ir [(3kT/b ) M u 

+ b<rf>(v - O X u) + (kT/b ) V In i//] (15) 

which yields 

M = ^[3«V'Mu> 

+ 6„ 2 (o T 0(v - ft x u )> + kT{u v V In i/>>]. (16) 

(1 + N^//)- 

The type of laminar flow shows up in v and fl. 

The bilinear coordinate averages in eq (16) can be derived 
from the diffusion equation, eq (11), by multiplication by f p 2 , 
^pVpi ' ' ' an( l integration over the whole space. 

5. Flow With Transverse Gradient 



[cr] = lim 



VI 



(r T n 



/i^\ In the case of laminar flow with transverse gradient v -h 

y(y,o,o), one obtains the set of linear differential equations 



Up 2 ) ~ Vp/Sk, - 2r p y(l + v p <p p /2f)(£ p r) p ) = 

(^pVp) - T p y(i> p <p p /2f)(i; p 2 ) - VyO + Vp<p P /2f)(Vp 2 ) = 

(Vp 2 ) ~ Vp/3\ p + r p y(v p <p p /f){t; p r) p ) = 

(VpCp) + Tpy(v P <pp/2f)(tp£;p) = 

OS) " *p/3A, = 

(£*>&>> + T p y(l + VpVp/Zf) (t)p£ p ) = 



-r P d(€ P 2 )/di 
~Tp d<i P 7j p )/dt 
-r' p d( Vp 2 )/dt 
-t'p dfy)p£p)/dt 
~r p d(£ p 2 )/dt 
-t' p d(^ p )/dt 



(17) 



The steady state solutions which apply after the complete 
decay of transient phenomena read 



UpVp) = 



i 
l 



l + 



yo 2 Tp(T P +T p )/2 

1 + icorL 



1 + 



1 + 2ia>Tp 



V£p_ eUot 



3/Xp 1 4- i(x)Tp 

I 



{vi ).-L\,- y '' T ' iT - - T ' ,/2 ( 

3^ P 



1 4- io)T p 



?io)t 



1 + 2i(or p 



(18) 



(Cp 2 ) = — 
3/Xp 

(VpCp)= (Cp€p) = 
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The averages are cut off beyond the lowest power y which is 
needed later in zero gradient expressions for intrinsic viscos- 
ity, normal stress difference and birefringence. These aver- 
ages have to be inserted in the expression for the intrinsic 
stress tensor. 



[&] = (N/M)[3kT(u T tAu) 
+ (yb 2 /2)(u T (f>(r), £o)> - kT}] X (1 4- NO//)" 1 (19) 

yielding the frequency dependence of intrinsic viscosity 

[ij] = W. - HvTo 

= lim — = km - — = Lo"i2r/T7s7 



c^O CT7 S c^O CO" 12s 



/fr y fSfJLp 

M VsP = 



1 + ico(Tp - T p ) 



1 + 



/ 



/?r 



Mi7 s p 14- icoTp 



RT 



[I + IV - i(II - III)]/co 



tan 8 V = (II - III)/(I + IV). 



(20) 



A rather similar but not identical expression applies to 
intrinsic streaming birefringence 



[A/*L* - [An]'. - [A,]". 

An* - An s An* - An s 
= lim = lim 



c-+o cn s T7 s y c -^o cn s o- 12s 



M-q sP =i y I \ 



RT v 2(jl p 

fs P 

4tt (n 2 + 2\2/V ai - a 2 ) 






5 \ 3/i / Mr) s 1 4- icoTp 

Mr) s 



tan 8„ = II/I. 



(21) 



In the case of dynamic birefringence An is the difference 
between the refractive indices in the diagonal direction in the 
first and second quadrant. The extinction angle X = 45 as 
long as the gradient amplitude y is small enough. The 
birefringence [An]JVIir) s /KRT and the phase angle 8 n are 
plotted in figure 3 for the same values of Z, if If and h* as in 
figure 2. 



The absolute value of intrinsic viscosity [r]]J[frj s /RT and the 
phase angle 6^ as function of cor l are plotted in figure 2 for Z 
= 100, <p/f '— 2, h* = .4 according to old and new theory. 





FIGURE 2. Relative intrinsic viscosity [t)]J[?)]o and phase 
angle d v for Z = 100, h* = .4, and (f/i = 2 as function of 
(OT l according to new (full line) and old (broken line) theory. 



1 

log (cut,) 

FIGURE 3. Relative intrinsic streaming birefringence 
[AnJ^/fAnJo intrinsic accoustic birefringence [An] ac and phase 
angle d n for Z = 100, h* = .4, and ip/f = 2 as function of 
0)T 1 according to new (full line) and old (broken line) theory. 
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The meaning of the symbols in eqs (17) to (21) is as 
follows: 



b 2 /6D k p 



t v ~ r v 



K 



T p (l + v P (p P /f) 

T P Vp<Pp/f 



4tt (n 2 + 2Wa 



5 \ 3n ) \ kT 
I = 2wr p /(1 + co 2 r; 2 ) 

ii = ^wv; /(i + w 2 r; 2 ) 

III = ^0) 2 T p (t' p - T p )/(\ + GJ 2 t' 2 ) 

IV - 2 «> 3 T P T P (t p - r p )/(] + orr^) 




300 



Figure 4. Second Newtonian intrinsic viscosity [rjjoc as 
function of number of links Z and hydrodynamic interaction h 
as parameter for {p\i — 2. 

New (full line) and <>l<l ll.mk.-n line) theory. 



The sums III and IV go to zero for vanishing internal 
viscosity <p — > 0. In this limiting ease the intrinsic viscosity 
and streaming birefringence are proportional to each other. 
Here/? is the gas constant, K is the rheooptical coefficient, // 
is refractive index, 8 is the phase angle between the flow 
gradient and viscosity or birefringence, OL\ and a 2 are the 
optical polarizabilities of the link in the directions parallel 
and perpendicular to the link respectively. 

Both definitions of intrinsic viscosity and birefringence in 
eqs (20) and (21) tend to make the intrinsic values as much as 
possible independent of the specific properties of the solvent, 
i.e., of 7] s and n s which influence the orientational forces, 
rigidity and optical anistropv of the dissolved macromolecule 
[52]. 

It is important to note that as a consequence of the new 
terms III and IV, i.e., (i7 P 2 ), caused by the introduction of 
internal viscosity the expressions for intrinsic viscosity and 
streaming birefringence, eq. (20) and (21), respectively, are 
not proportional to each other. Hence the rheooptical law 
does not apply to such a model. In particular [Aft]^ cannot be 
written as K[r)]co ~ [cr] 12 . 

The proportionality, however, still holds in the first New- 
tonian region, (x) —> 0, where all the terms but I/a) go to zero. 
This means that the intrinsic viscosity [17] and the Maxwell 
constant [A/z] and hence the stress and the optical tensor are 
proportional to each other for y — > and a) —> 0. This is the 
case in most applications of flow stress mapping by means of 
streaming birefringence of dilute polymer solutions. 

With increasing frequency and finite internal viscosity, 
[An] w goes to zero while [tj]^ tends to finite value. The 
rheooptical law breaks down completely. In this second 
Newtonian range the intrinsic viscosity is independent of 
frequency. One derives from eq (20) 



h]« = T7~ 2 



The limiting values [vUMirjJRT forZ between 3 and 300 
and h* = . 1, .2, .3, .4 are plotted in figure 4 for the new (full 
line) and old (broken line) theory. One sees that the absolute 
values for v p = 1 and V p l) = k pZ /k pR differ by less than a 
factor of 2. But their dependence on Z, i.e., on molecular 
weight, is just the opposite for h* = .4. In the incorrect 
formulation of the old theory [7]]^ slightly increases with Z 
while in the correct formulation of the new theory it decreases 
withZ. Such a dependence on Z is much more in agreement 
with experimental data on polystyrene in Aroclor [44]. The 
second Newtonian viscosity [tjJoc is 14.3 cm 3 /g if one goes 
with M from 20,000 to 860,000. The old theory yields a 
steady increase of [r)] x with M in sharp disagreement with 
these data. 

The often used real and imaginary part of intrinsic shear 
modulus 



[C]* = [G]' + i[Gj = fcwfcfaE 



= t[ct 12 ]* = lim 



cy 



(24) 



Mt7 sP =i 1 +f/v p <p p 



(23) 



are plotted in figure 5. The difference between the old and 
new theory is relatively small forZ = 100 but would be larger 
forZ = 300. 

The consequence of the non- vanishing second Newtonian 
viscosity, [rjjoc =£ 0, is the linear increase of [G]' with o>t, in 
the assymptotic high frequency region. Such a behavior is in 
perfect agreement with experimental data on polystyrene in 
highly viscous Aroclor [44-51]. These data together with the 
zero gradient intrinsic orientation data of streaming birefrin- 
gence of polystyrene in solvents of increasing viscosity [53) 
constitute the main support of the theory. 
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A* - Ae~^ 



Figure 5. Intrinsic shear moduli [G]' and [G]" for Z = 
100, h* = .4, and <p/{ = 2 as functions of cojj according to 
new (full line) and old (broken line) theory. 



RT 



y P 2 



- + 



M /P £i 1 + icor; \1 + 2ia)T p 

[1 + 2«*t£ (2r' p - Tp )] 



«T . 2 y 2 - iaifo + 2t p + 4co 2 t; 2 (t; - r p )] 
M 7o PI Tp " (1 + o> 2 t; 2 )(1 + 4co 2 t^ 2 ) 



(26) 



is proportional to the square of the amplitude of the oscillat- 
ing gradient as in the case of no internal viscosity. The value 
(A/y 2 )(M/RT), the phase angle 8 U and the frequency inde- 
pendent displacement 



(B ■ cos 8 2 /y z )(M/RT) = 2 



~ 1 



(27) 



are plotted in figure 6 as functions of (or 1 for Z = 100 and h* 
= .4. Both quantities A and B are going to zero with 
increasing frequency. In contrast to intrinsic viscosity the 
limiting first normal stress difference, at (O — » oc , does not 
become finite by introduction of internal viscosity although it 
goes to zero more slowly, as (O l instead of as co -3 . 



The intrinsic first normal stress difference turns out to be 



[0-11-0-22]* 



rt -<~ ■«*>«*> 



hi I / 



(1 + Vp(p p /f) 



*L 2 y ^ 

M 7o £i 1 + ia>r p 



(25) 



1 + 



1 + 2ia)Tp 



EAJtot 



= A*e 2ioit + B* 

while the second normal stress difference vanishes in all 
isolated necklace models with ideal elastic links independent 
of internal viscosity. In applying eq (25) one must not forget 
that the phase angle 8j of A* is dependent only on the factor 
of exp (2ia)t) and that the angle § 2 only reduces the constant 
vertical displacement to B cos 8 2 but does not yield any phase 
shift. 

The term 1 in the parenthesis, i.e., /?*, keeps the first 
normal stress difference positive during the whole period up 
to very high frequencies. The oscillation is taking place with 
twice the frequency of the flow field. All these effects are very 
much the same as in the case of no internal viscosity. The 
difference is mainly in the replacement of T p by 7 p and the 
term proportional to r' p 
oscillating term 



r p . The amplitude A of the 




FIGURE 6. The coefficient B, the sum B + A, the differ- 
ence B — A of the relative intrinsic normal stress difference 
[°-n-0-2 2 L/[o-n-cr 22 ]o = A ei 2i *"- i81 + Bfor Z = 100, h* 
= . 4 and <p/f = 2 as functions of (*)T 1 according to new (full 
line) and old (broken line) theory. 
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6. Flow With Longitudinal Gradient 

y(— jc/2, — y/2,z) the intrinsic stress 



In jet flow, v 
tensor 



3RTy 



T p (l - (OTp 

1 + a>V 



+ 



T» Tv 



M 



1/ 



[:U-7'<uMo> 






(I/oj + E (r; - T„)/3 - ;il/a>). 



+ {yb 2 /2){» </>(-£, -»}, +2t,)> + ATI] 
x (1 +N <£//)-' 



(28) 



The additional term 2j ( t /> — t p)/3 * s independent of O) and 
hence represents the Trouton viscosity in the second Newton- 
ian range where I/o) and II/co disappear. As in the case of 
is independent of coil rotation in flow because the flow field conventional intrinsic viscosity, eq (20), the finite value at a> 



has no rotational component, ft = 0. The bilinear coordinate 
averages are derivable from the set of differential equations 

(1 + Jt p ){^) - v p /3\ p = -T p d(&)/dt 

(i + jTvKipVp) = - T v <i(€pVp)/<Ii 

(1 + YTpKVp 2 ) - fp/'^K = -r'p <1(Vp 2 )/'Ii 

(29) 

(1 - yr v /2){j] l> Q = -t' v <l(vp(,p)/<li 

(i - nMp 2 ) - v P /$K = -t'„ d{tf)/dt 
(] - yT p /2X{„&> = -t'„ d(Up)/dt 

Under consideration of the svmmetry of flow field one 
obtains 

1 

3fi p 



<fp 2 > = (rip 2 ) = — exp {-[/ + (yoT p /iio)e*"]/T p } 



Jo 



sxp|>+ (yoTjiw'p^'^dx 






f 

Jo 



exp [x - 2{y r p /io)r p )e UtiT '^ dx 



(ZpVp) = (VpUp) = (Zp^p) = °- 



(30) 



<^ 2 > 
(ip 2 ) 



I - JoTp 



1 4- ia)T p 



3fi p 

1 / . e**' 

— ] + 2y T p -— — 

3/Ltp \ 1 + iWT p 



(31) 



The intrinsic Trouton viscosity turns out to be 

ML 



(32) 



[cr 33 - <r n ]*/yri 8 

/^T? s P=l 



00 is a consequence of internal viscosity, i.e., of partial 
coil rigidity. The frequency dependence of [77]^ and phase 
angle 8 iv are plotted in figure 8 for Z = 100, h* — .4, and 

The intrinsic birefringence reads 



[Ant = K-^r 2 :v P <&> 2 - & 2 > 



(33) 



which for small amplitude reduces to 

r ,* 3A/?r ^ , , x 3K/?r 

[A/i]L = -— • Z t,/(] + icori ) = (I - ill) 



tan d bl = 11/ I 



(34) 



These averages still contain the transient which, in the 
general case, cannot be easily separated from the stationary 
solution reached after the transient has tapered off (fig. 7). 

The separation can be performed if y T p /(or' p is so small 
that one can replace the exponential function by its linear y q 
expansion. In such a case one obtains for the stationary 
solution 



This expression differs from the streaming birefringence in an 
oscillating flow field eq (21), only in the factor 3. The non- 
proportionality between viscosity and birefringence is again 
the consequency of the additional internal viscosity term in 
eq (33). 



2.0 



1.5 - 




t/277-T'p 

FIGURE 7. Time dependence of (£ p ~) = (i7 P 2 ) according 
to eq (30) plotted versus i/27TT r p for the special case ipr' p = 1 
and a = y T p/ WT 'p showing the short duration of transient 
effects and the rapid approach to the asymptotic periodic 
solution. 

The stabilization is slower for low o>t' p . 
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Z=100 
h # =.4 
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-2-101234 
log tUT| 

FIGURE 8. Relative intrinsic Trouton viscosity [v)]ij\v)\m 
and the phase angle 3^ for Z = 100, h* = .4, and <pj{ = 2 as 
function of ojt^ according to new (full line) and old (broken 
line) theory. 



Acoustic birefringence is in many respects closely related 
to birefringence in an oscillating jet flow. The main differ- 
ence is the absence of lateral contraction as represented by 
— yxl'2 and —yy/2. In contrast to jet flow with constant 
specific volume (incompressible liquid) the volume element 
subjected to an acoustic wave is periodically compressed and 
expanded. That means a constant (£ p 2 ) yielding in eq (34) the 
replacement of the factor 3 by 2. 

One has 

r i An 

LArcJ*. = hm — — - 



N 



0)T 



m ut (i + oj 2 t 2 ) 1 

z 



sin [(x){t - zq/co) - 8 ac ] 



N 
M 



N 






p=l 1 + idJTj 



= - K ac (I - HYSe^-** 10 ** 
M 



The difference between acoustic birefringence according to 
the correct new and the incorrect old values v p can be seen in 
figure 3 where [An] ac is plotted versus cor l forZ = 100, h* = 
.4 and <p[f = 2. The phase angle 8 ac is identical with 8 n and 
Sin- 

7. Conclusions 

The paper presents the calculation of most of the intrinsic 
rheological and rheooptical effects of linear homopolymers in 
an oscillating flow field which may be explored experimen- 
tally. In the case of rheooptical effects only the intrinsic 
birefringence of the polymer is included. The form-birefrin- 
gence is completely neglected. The same applies to the 
influence of large side groups which may effect independ- 
ently the optical anisotropy of the segment and its frequency 
dependence. 

The introduction of the appropriate V p — 1 values instead 
of the old values v p — A pZ /A. p # does not change drastically the 
effects depending on internal viscosity. As a rule the ratio 
between the new and old values of intrinsic viscosity, bire- 
fringence, and first normal stress difference is less than 2, at 
least in the range ofZ between 3 and 300. With higher Z the 
differences increase as a linear function of log Z. 

As already mentioned, the smallness of the difference is a 
consequence of the peculiar dependence of old v p(X) on p: 
much larger than 1 at small/? and smaller than 1 at high/). 
Hence the larger contributions in the former part of the sums 
are partially compensated by the smaller contributions in the 
latter part. 

The most important changes occur in the second Newton- 
ian intrinsic viscosity which is the most conspicuous conse- 
quence of internal viscosity. Here the dependence of [?)]x on 
Z is much less according to the new theory than it was in the 
case of the old theory. 

No calculation of the gradient effects was attempted be- 
cause one knows that in a flow with a finite velocity gradient 
the randomly coiled macromoleeule is deformed with a con- 
sequent change of all interbead distances. This yields a 
change of interaction tensor H which leads to a modification 
of all eigenvalues k p . The rest values X p used in this paper 
are only applicable to effects where an extrapolation to zero 
gradient is straightforward. This is the case with dynamic 
effects where one uses very small gradients and concentrates 
on the frequency dependence of the effects measured. The 
situation, however, is basically different in the non-linear 
range of gradient dependence of excess stress and optical 
tensor. 



tan § a 



H/I 

47T In 1 + 2\2/ 2 \i/2 > 



5 \ on 
l ac = pB 2 c ac /2. 



Pc*a 



(«l - oc 2 ) 



(35) 



Here z is the location of the center of hydrodynamic resist- 
ance of the macromoleeule, B is amplitude, l ac is intensity 
and c ac is propagation velocity of acoustic wave, and p is the 
density of solution. The frequency dependence of acoustic 
birefringence and phase angle is exactly the same as in the 
case of oscillating jet flow. 
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